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ABSTRACT 

GJ 758 B is a cold (^600K) companion to a Sun-like star at 29 AU projected separation, which was 
recently detected with high-contrast imaging. Here we present photometry of the companion in seven 
photometric bands from Subaru/HiCIAO, Gemini/NIRI and Keck/NIRC2, providing a rich sampling of 
the spectral energy distribution in the 1-5 /im wavelength range. A clear detection at 1.58/im combined 
with an upper limit at 1.69/im shows methane absorption in the atmosphere of the companion. The mass 
of the companion remains uncertain, but an updated age estimate indicates that the most likely mass range 
is ~30-40 Mj U p. In addition, we present an updated astrometric analysis that imposes tighter constraints 
on GJ 758 B's orbit and identifies the proposed second candidate companion, "GJ 758 C", as a background 
star. 

Subject headings: planetary systems — brown dwarfs — techniques: high angular resolution 
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1. Introduction 

In recent years, a number of high-contrast compan- 
ions have been thermally imaged around nearb y stars 



(e.g. iMarois et al.l I2008t lLagrange et al 



2010). 



One 



interesting companion is GJ 758 B (Thal mann et al 
20091 hereafter PI). Its combination of a Sun-like 



parent star (spectral type G8V) only 15.5 pc away 
(IPerrvman et al.ll 1997b . the close proximity to the star 
(projected separation 29 AU), and low surface temper- 
ature (^600 K) make it one of the most "planet-like" 
objects available for direct study, acting as a laboratory 
for current planet formation and evolution theories. 

PI presented two epochs of i/-band imaging of 
GJ 758 B that allowed for proof of common proper 
motion and a first estimation of physical and orbital 
properties. Furthermore, a candidate second com- 
panion, tentatively called "GJ 758 C", was found in 
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one epoch. Recently, ICurrie et al. (l2010l) published 
L'-band data that confirmed the estimated tempera- 
ture range of B but did not detect "C". Both publi- 
cations deduced similar results for best-fit orbital pa- 
rameters (a ~ 50 AU, e ~ 0.7). Both also noted 
that the mass range is ^10-40 Mj Up if the full range 
of main-sequence ages for a Sun-like star is consid- 
ered. In this work, we present photometric coverage of 
GJ 758 B in the near infrared ( JHK C L'M bands) with 
additional photometry in the methane-sensitive CH4S 
and CH4L narrow-band filters (see Table Q] for fil- 
ter specifications), collected with high-contrast imag- 
ing techniques at Subaru/HiCIAO, Gemini/NIRI, and 
Keck/NIRC2. This is used for updating the estimations 
of physical and orbital parameters of GJ 758 B. We re- 
port the re-detection of "GJ 758 C" and the conclusive 
identification of it as a background star by proper mo- 
tion. 

2. Observations and Data Reduction 

i/-band observations of GJ 758 B were obtained 
with Suba ru on November 1 1, 2009. making use 
of AO 188 dHavano et alJ EoToh coupled with HiCIAO 
(iHodapp et al l 120081) as part of the SEEDS survey 
dTamural2009l) . The field of view (FOV) was 20" x 20" 
with a 07010 pixel scale. The total integration time 
was 906 s. 

Observations in J-, H-, CH4S-, CH4L-, and 
if c -band were obtai ned at Gemini Nor th with the 
Altair AO system dHerriot et all l2000h and NIR1 



( Hodapp et al] 2003 ). The observations were carried 
out on 5 different nights between April 27 and May 
08, 2010, with a 0'.'022 pixel scale over a 22"x22" 
FOV. The total integration times were 900 s, 1800 s, 
3810 s, 2220 s, and 2400 s, for the J-, H-, CH4S-, 
CH4L-, and K c - band, respectively. The nearby star 
HD 226294 was observed before or after each ADI 
dataset in the same filter to be used as a photometric 
reference. 

We imaged GJ 758 in the V- and Af s -bands 
on August 6, 2010 using the Keck II AO system 
(Ivan Dam et all |2004 and NIRC2. The 0"300 di- 
ameter coronagraphic mask was used for the V ob- 
servations to prevent saturation. Images without the 
mask were also taken, to calibrate the brightness of the 
star. For all observations, we avoided the bad quadrant 
in the NIRC2 detector and used the 768x776 subar- 
ray mode with a O'.'OIO plate scale. The integration 
times were 1800 s and 1995 s for V and M s , respec- 



tively.Three photometric standard stars (HD 162208, 
HD 161903, Gl 748) were observed in both bands. 

All data wer e taken using angula r differential imag- 
ing (ADI, e.g. iMarois et al. 2006) with typical field 
rotations of ^30°, a nd the images were red uced with 
the LOCI procedure dLafreniere et alj|2007l) . using an 
IDL implementation adapted for each of the instru- 
ments. For the HiCIAO data, the same procedure as 
described in PI was used, including correction for the 
partial subtraction of companion flux during the LOCI 
reduction. The only practical differences for NIRI and 
NIRC2 were the image registration procedures, where 
the NIRI registration was performed on the basis of 
cross-correlation of the diffraction spiders, and NIRC2 
registration was caclulated by centroiding on the stel- 
lar PSF core, which was non-saturated in the M-band 
data and non-saturated behind the semi-transparent 
mask in the L'-band data. Also, for the NIRC2 data, 
high-pass filtering was applied at the same time as the 
subtraction of a radial profile from the stellar PSF in all 
images, to remove low-frequency spatial variations in 
the hi gh thermal backg round of the L'- and M s -band 
(e.g. Janson et alJl2008 ). Fluxes of all observed targets 
and standard stars were extracted using aperture pho- 
tometry. The reduced images are shown in Fig. Q]and 

Fig.m 

3. Photometric Analysis 
3.1. Calibration 

Photometric calibration is of central importance for 
these data, due to the fact that the JHK S -photometry 
of GJ 758 A has been flagged as unreliable in 2MASS 
( Skrutskie et al.ll2006h . and furthermore, for the broad- 
band filters it is always saturated in NIRI images, even 
for the shortest available integration times in subarray 
mode. Hence, we used the standard star for photo- 
metric calibration in all filters, with the exception of 
CH4S. Those data were taken on Apr 29 under photo- 
metrically unstable conditions, so that the standard star 
photometry couldn't be trusted. Instead, we calibrated 
the CH4S flux based on the fluxes of the background 
stars in the GJ 758 field in CH4L and CH4S, under the 
assumption that they are equally bright in both filters, 
which is reasonable since they are stars and cannot ex- 
hibit methane absorption. Four background stars were 
used with a dispersion of 0.19 mag. The brightness of 
the background stars in CH4L was calibrated against 
the standard star. The resulting absolute magnitudes 
of GJ 758 B are summarized in Table [T] 
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Fig. 1. — High-contrast imaging of GJ 758 B with ADI/LOCI. In all panels, the star is located approximately in 
the upper left corner, North is up, and East is left. The location of GJ 758 B is marked with a solid arrow, that of 
"GJ 758 C" with a dashed arrow, (a) HiCIAO ff-band image (Nov 2009). GJ 758 B is clearly visible, whereas "C" is 
marginally detectable at ^3cr. (b-f) Gemini/NIRI images in five filters (Apr/May 2010). (g-h) NIRC2 images in L'- 
and M s -band (Aug 2010). 
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Fig. 2. — Combination of Gemini images, (a) Co-add of J, H, K c , and CH4S, after normalizing the flux levels such 
that GJ 758 B appears approximately equally bright in each image. The "C" source is detectable, (b) A map of the 
pixels that exhibit a S/N ratio of >+l in all four images used for (a). Only GJ 758 B and "C" pass the test, indicating 
that they are both real on-sky objects. 
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GJ 758 A has no published L'-band or Af s -band 
magnitudes, thus we use the standard stars to calibrate 
those values. All three standard stars give consistent 
results within 0.1 mag, yielding Mry = 3.7 ± 0.1 
mag and Mms = 3.6 ± 0.1 mag. During the sci- 
ence exposures, GJ 758 A was unsaturated in the Ms- 
band and so could be used to evaluate the brightness 
of GJ 758 B. In the case of L', the star was behind 
the semi-transparent mask, which makes it unreliable 
for photometry, hence for this case, standard star pho- 
tometry was used directly to calibrate the brightness 
of GJ 758 B. The results are in TableQ] There are two 
cases where special circumstances apply: For CH4L, 
the flux at the position of the companion is visually 
unconvincing as a point source, and more reminiscent 
of residual noise from the stellar PSF. It is therefore 
possible that GJ 758 B is systematically overestimated 
at this wavelength, such that the quoted CH4L flux 
should rather be regarded as an upper limit rather than 
the true brightness. In the case of M s -band, there is 
no visually clear signature of the companion, but aper- 
ture photometry was performed on the same location 
as the L'-band detection. This yielded an excess flux 
at the ^2<t level, hence this is a plausible flux of the 
companion, since its position is known a priori. Aside 
from CH4L and M s , all detections are ^10er confi- 
dence. For NIRI, the photometric error that we cal- 
culate is dominated by an uncertainty in the linearity 
behaviour of the detector at low counts/short integra- 
tions, which we estimate could impose an error of at 
most 20%. 

3.2. Interpretation 

GJ 758 B exhibits clear methane absorption. This 
is generally expected for objects in the late T-type 
range, although to a decreased extent if the atmosphere 



is in chemical non-equi librium (e.g. Burg asser et al 



2006; Fo rtnev et al.ll2008l) . GJ 758 A has a metal 



licity of approximately +0.2 dex (e.g. Holmber g et al 



2009; iKospal et al . 2009). Assuming that the compan- 
ion has the same composition, this super-solar metal- 
licity is likely at least partly responsible for the high 
flux in K c . Since t he commonly used COND models 
dAllard et alJbOOlHBaraffe et al.ll2003l) do not contain 
any non-solar abundances at the relevant temperature 
range, we do not use them for model comparison. In- 
stead, we use models based on iBurrows et al. (2006) 
extended to colder temperatures (Hubeny & Burrows 
in prep.). With our new adopted age range discussed 
below, the feasible values of log g for the companion 



are all close to 5.0, hence we use this as a baseline 
value in our fitting. The comparison was done by fit- 
ting model spectra to the flux densities of the photo- 
metric points, which were derived from the magnitude 
values using the NIRI and NIRC2 filter transmission 
curves. All fluxes were normalized to a distance of 
10 pc. The best-fit value was acquired using the mini- 
mized x 2 m me photometric bands as the quality met- 
ric. The models cover temperatures of 500-700 K in 
steps of 50 K, log g of 4.0-5.0 in steps of 0.5 dex, and 
metallicities of +0.0 and +0.5 dex (note that the full 
grid is not covered - in particular, no super-solar abun- 
dace models exist below 600 K). The best-fit temper- 
ature of this procedure is 600 K, in agreement with 
previous results. A metallicity of +0.5 dex fits better 
than +0.0 dex, which is consistent with the metal en- 
richment expected in the system. However, since sur- 
face gravity and metallicity are largely degenerate for 
the wavelengths covered by our study and the models 
are very uncertain, we do not attempt to optimize the 
fit with respect to either of these quantities. Indeed, 
the models are known to be unable to reproduce vari- 
ous features of other c ool companions to stars (e.g. the 
HR 8799 planets see iMarois et al.ll2Q08t Ijanson et al 



2010; lHinz etal.lE010l:lBowleret al.ll2010l) . so deriv- 



ing absolute values from model comparison is of lim- 
ited relevance. Four examples of model fits are shown 
in Fig. [3] 

We also make a comparison with field brown 



dwarfs fromLeggett et al. (2010). For this purpose, we 
use the J-band absolute brightness M j = 17.58 ± 0.20 
mag, and the color J — H = 0.58 ± 0.28, mag, since 
this should be relatively insensitive to metallicity and 
gravity. The compa nion fits very well to the tempera- 
ture sequence of the lLeggett et al. ( 2010b field dwars in 
an HR-diagram, among the latest-type brown dwarfs 
(T8-T9). The J — H color alone with uncertainties 
places the companion firmly beyond T5. A T8-T9 
spectral type fits well to the temperature of ~600 K 
derived from model comparison. 

In PI, we adopted an age range of 0.7-8.7 Gyr, on 
the basis of a variety of age indicat ors. The 0.7 Gyr ag e 
was based on isochronal fitting bv lTakeda et al.l (120071) . 
However, other isochronal fits have been performed 
yielding substantially higher ages of several G yr (e.g. 
Valenti & Fisched 120051; llolmberg et al.ll2009h . This 



indicates that isochronal fitting is inadequate for dating 
GJ 758, hence we revise our age estimates by remov- 
ing the isochrone method, which yields a residual age 
range of ^5-9 Gyr (keeping estimates based on activ- 
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ity and rotation). As an additional age determination 
test, we have analyzed a high-resolution (i?~3 1,500) 
spectrum of GJ 758 from the Apache Point Observa- 
tory 3.5m. We detect no clear evidence of the Li I 
6708 A doublet at an upper limit of ^4 mA equiv- 
alent width. If Li had been detected, it would have 
indicated a young age, but the non-detection does not 
provide a strong lower limit on the age, as it only sug- 
gests that GJ 758 sh ould be older than the ~ 100 Myr 
old Pleiades cluster dMaldonado et al. 2010l) . We re- 
iterate that age determination of main-sequence Sun- 
like stars is highly uncertain, and that the adopted age 
range of ^5-9 Gyr should not be considered definitive, 
but merely represents the range of mean ages from 
methods that are considered sufficiently reliable. The 
corresponding range of companion masses from evolu- 
tionary models is ^30-40 Mj up . The most promising 
avenue for getting a reliable age estimate of GJ 758, 
and thus a better constraint on the mass of its compan- 
ion, is likely asteroseismology. 

4. Astrometric Analysis 

4.1. Data points 

In PI, we demonstrated the common proper motion 
of GJ 758 and GJ 758 B on the basis of two epochs 
of observation, May 3, 2009 (El) and August 6, 2009 
(E2). In this work, we include three additional epochs 
in our analysis: Our data from November 1, 2009 
(E3) and April 29, 20 10 (E4), as well as the May 27, 
2010 data point from ICurrie et alj d2010b (E5). The 
Keck/NIRC2 data from August 2010 (E6) have a nar- 
row field of view excluding the 5-7 known background 
stars that were used to fine-tune the pixel scales and ro- 
tation angles in our previous observation. Thus, they 
do not deliver sufficient astrometric accuracy. 

Figure^ shows the resulting positions of all nearby 
point sources relative to GJ 758 at E1-E4, as well 
as the data point for GJ 758 B in E5. The E4 data 
points represent the mean positions of the sources in 
the four datasets that reveal GJ 758 B, i.e. the J, 
H, K c , and CH4S data. The error bars for the Hi- 
CIAO data (first three epochs) are based on the pixel 
size of the HiCIAO camera, 9.5 mas. For the Gem- 
ini data, the pixel size is 22 mas, yielding an error of 
22 mas/ V^- = 11 mas for the combined data points. 
These errors are consistent with the scatter of the back- 
ground star data points around the projected motion 
path, given the expected contributions from the proper 
motions of the individual background stars. For the 



E5 data point, we assume an isotropic error of 10 mas 
to repres ent the anisotropic e rror bars of 5 and 15 mas 
shown in lCurrie et alJ (120101) . 

We observe that GJ 758 B pursues a trajectory to 
the northwest consistent with orbital motion relative to 
GJ 758, clearly setting it apart from the background 
star trajectory, which is dominated by GJ 758's known 
parallactic and proper motion. The source tentatively 
referred to as "GJ 758 C" in PI is found to follow 
the background star track; while this was still indis- 
tinguishable from orbital motion at the 2a level in 
November 2009, the May 2010 data are unambiguous. 
This also allows u s to identify the candidate signal in 
Curri e et al. (2010) as spurious 



4.2. Orbital Monte Carlo simulation 

With a total of A = 10 scalar parameters (2 coordi- 
nates x 5 epochs), the astrometric data on GJ 758 B 
is now extensive enough to fit synthetic orbital so- 
lutions generated by a Monte Carlo simulation with 
the least-squares method rather than with the simpli 



fled a pproach previously used in PI and ICurrie et al. 
(2010h. However, since the curvature of the orbit is 



not yet measurable, the benefit of this improvement is 
limited. 

We generate a large number (>10 6 ) of orbital tra- 
jectories with random values for eccentricity e, in- 
clination i, argument of periastron uj, and longitude 
of the ascending node £1 The distributions are pre- 
sumed to be flat, except for the inclination, where 
larger angles are favored proportionately to sin? in 
order to represent their higher geometric likelihood. 
The two remaining orbital parameters, the semimajor 
axis a and the mean anomaly at epoch Mq, are im- 
plicitly chosen by defining an anchor point (xa,Ua) 
in the projected image plane where the companion is 
located a given epoch £a- In order to achieve a high 
production rate of valid orbital solutions, we choose 
tA to be the mean of the five observational epochs, 
(iobs), and generate (xa,Ua) randomly in a box of 
20 mas x 20 mas centered on the mean astrometric co- 
ordinates ((a; b S ), (y bs)). 

In order to evaluate an orbital solution for consis- 
tency with the data, we determine the predicted po- 
sition of the companion at the five epochs and calcu- 
late the x 2 deviation. The minimum best fit achieves 
X 2 = 2.09. We select the set of orbits with x 2 < 3.0 w 
■f 1 to represent the "good fits". 
Figure |4p illustrates the distribution of semima- 



Y 2 . 
A mm 
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Wavelength (/i.m) Wavelength (/xm) 



Fig. 3. — Photometric analysis of GJ 758 B. The red plus signs show the measured flux values normalized to a distance 
of 10 pc for the seven filter bands, whose wavelength domains are marked with black bars. The orange curves represent 
model spectra for different assumptions of T e s , log g, and metallicity, and the black diamonds are the resulting flux 
levels in the filter bands. 
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jor axis a, eccentricity e and inclination i for 2516 
"good fit" orbits. The solutions appear to lie in a two- 
dimensional manifold in the three-dimensional param- 
eter space. We note that more stringent fitting require- 
ments (e.g. x 2 < 2.5) do not relevantly reduce the 
size of this manifold, showing that this spread is due 
to model degeneracy rather than fitting errors. The 
degeneracy represents the fact that only the projected 
position of GJ 758 B can be tracked, leaving the line- 
of-sight component of its position and velocity unde- 
termined. In order to break this degeneracy, it is nec- 
essary to measure the curvature of the projected orbit. 
Given the typical predicted orbital periods of several 
centuries, this requires monitoring on time scales of at 
least a decade. 

The numerical results are summarized in Table Q] 
We note that the w eighted median ecc entricity has 
dropped from 0.73 in lCurrie et alJ (1201 0) to 0.56; this 
is because the new epochs E3 and E4 suggest that the 
previous epochs were overestimating the orbital veloc- 
ity. 

We thank David Lafreniere for providing us with 
the source code for LOCI, Eric Mamajek for use- 
ful discussion, the staff at Subaru, Gemini and Keck 
for their support, and an anonymous referee for use- 
ful suggestions. The Gemini time was allocated by 
NOAO. The Keck telescope was funded by the W.M. 
Keck foundation, and time was allocated by NASA 
through partnership with Caltech and UC. We ac- 
knowledge the cultural significance that the summit of 
Mauna Kea has to the indigenous Hawaiian commu- 
nity. Part of the research was supported by NSF grants 
AST 1009203, AST 0802230, AST 1009314, the AAS 
Chretien grant, a Grant-in- Aid for Specially Promoted 
Research, the Mitsubishi Foundation, and JPL, Caltech 
under NASA contract. 

Facilities: Subaru (HiCIAO, A0188), Gemini 
North (NIRI), Keck (NIRC2). 



Table 1 : Numerical Results 



Photometry (mag) 


GJ 758 A 


GJ 758 B 


Mj (1.15-1.33 ,um) 
M H ( 1.49- 1.78 Aim) 
M Kc (2.08-2.11 fjm) 
M L - (3.43-4.13 /xm) 
M M s (4.55-4.79 ,uni) 
M CH 4S (1.53-1.63 urn) 
M C H4L (1.64-1.74 Atm) 


3.7±0.1 
3.6±0.1 


17.58±0.20 
18.16±0.20 
17.12±0.20 
15.0±0.1 

> 10 Q + 0.5 

17.74±0.20 
>18.86±0.20 



GJ 758 B orbital parameters 


Weighted median 


68% interval 


Semimajor axis a (AU) 


44.8 


30.7-89.2 


Eccentricity e 


0.564 


0.416-0.712 


Inclination i (deg) 


43.6 


26.2-56.7 


Period P (yr) 


299 


170-843 



GJ 758 B's position relative to GJ 758 
Epoch ARA (") ADec (") Error (") Ref. 



El: May 3,2009 


-0.574 


-1.789 


0.010 


(1) 


E2: Aug 8, 2009 


-0.579 


-1.765 


0.010 


(1) 


E3: Nov 1,2009 


-0.597 


-1.751 


0.010 


(2) 


E4: Apr 29, 2010 


-0.609 


-1.735 


0.011 


(2) 


E5: May 27, 2010 


-0.616 


-1.716 


0.010 


(3) 



E6: Aug 6, 2010 — insufficien t astrometry — 

References. (1) PL (2) this work. (3) ICurrie et al.1 
( 2010h . Note that the error bars for astrometry in El 
and E2 were mistakenly listed as 5 mas in Table 1 in 
(1); however, the text and plots use the correct value 
of 9.5 mas. 
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Motion relative to GJ 758 in RA (mas) Semimojor axis a (AU) 



Fig. 4. — Astrometric analysis, (a) Motions of point-sources near GJ 758 across five epochs (E1-E5), measured 
relative to GJ 758's position. GJ 758 B exhibits common proper motion with its parent star as well as systematic 
orbital motion towards the northwest, whereas all other point-sources follow the expected trajectory for background 
stars (solid arrows). The object referred to as "GJ 758 C" in PI is unambiguously identified as a background star 
(motion highlighted by dashed blue arrows). The grey plus signs are la error bars. The circle marked as "PSF" shows 
the size of the resolution element in i//-band on HiCIAO. (b) Plot of eccentricity e against semimajor axis a for 2516 
orbital solutions with \ 2 < 3.0 generated by the Monte Carlo simulation. The orbit selection is biased according to 
the statistical weight (i>)/i>ob s as in PI. The inclination i is shown by color coding. The weighted median values of a 
and e are marked with a white plus sign. 
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